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Abstract 

Using the notion of the composita, we obtain a method of solving iterative func- 
tional equations of the form A 2 " (x) = F(x), where F(x) = ^2 n> o f{n)x n , /(l) ^ 0. 
We prove that if F{x) = Yl n >o f( n ) xn has integer coefficients, then the generating 
function A(x) = X^™>o a ( n ) xn i which is obtained from the iterative functional equation 
4A(A(x)) = F(4x), has integer coefficients. 

Key words: iterative functional equation, composition of generating functions, com- 
posita. 
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on : 1 Introduction 



For a given A(x) the n-th iteration is the function that is composed with itself n times. 
A°(x)=x, A l (x) = A(x), A 2 (x) = A(A 1 (x)), . . . , A n (x) = A(A n ~ l (x)) 



and denoted by A n (x). 

By an iterative functional equation we mean the equation, where the analytical form of 
a function A{x) is not known, but its composition with itself is known. 

Iterative functional equations arise in various fields such as fractal theory, computer sci- 
ence, dynamical systems, and maps. In particular, a challenging iterated functional equation 
A(A(x)) = x 2 — 2 circulate as a test of one's intelligence. 

Despite their prevalence, they are very difficult to solve, and few mathematical tools 
exist to analyse them. In the general case, solving the equations involve deep mathematical 
insight and experimentation with different substitutions and reformulations. 

In the literature you can find extensive investigations concerning iterates and functional 
equations, classical books [1, 2, 3, 4] and recent surveys and papers [5, 6, 7, 8, 9]. 

In this paper, we propose using the notion of the composita to identify and solve such 
iterative functional equations, where F(x) = ^2 n>0 f(n)x n , f(l) ^ 0. 
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2 Preliminary 



In the paper [10], the authors introduced the notion of the composita of a given ordinary 
generating function F(x) = J2 n >o f( n ) xn - 

Definition 1. The composita of the generating function F(x) = Y2 n >o f{ n ) xn is the function 
of two variables 

F A (n,k) = £ /(A!)/(A 2 ).../(A fe ), (1) 

where C n is a set of all compositions of an integer n, tt^ is the composition Y2i=i \ — n m to 
k parts exactly. 

The generating function of the composita of F{x) is equal to 

[F{x)f = Y,F\n,k)x n . (2) 

n>k 

The composita is the coefficients of the powers of an ordinary generating function 

F A (n,k) := [z n ]F(x) k . 
For the case k = n the composita F A (n, k) (see the formula (1) is 

F A (n,n) = f(ir. (3) 
In tabular form, the composita is presented as a triangle as follows 





/(I) 




/(2) 




/ 2 (i; 




^ 2 




r 4,2 


F A 



/(3) ^3 A 2 / 3 (1) 

m f a 2 f a ni) 

f{n) F A 2 ... ... d /"(l) 

Comtet[ll, p. 141] considered similar objects and identities for exponential generating 
functions, and called them potential polynomials. In this paper we consider the more general 
case of generating functions ordinary generating functions. For more information about the 
composita you can see [10, 12, 13]. 

In the following theorem we give the formula for the composita of a composition of 
generating function, which was proved by the authors in [10]. 

Theorem 2. Suppose we have the generating functions A(x) = ^2 n>0 a(n)x n , G(x) = 
S n >o 9i n ) xn ' an d their compositae A A (n,k), G A (n,k) respectively. Then for the compo- 
sition of generating functions F(x) = G (A(x)) the composita is equal to 

n 

F A (n, k) = J2 AA ( n > m)G A {m, k). (4) 

m=k 
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3 Iterative functional equation A(A(x)) = F(x 

The simplest case of an iterative functional equation for an unknown function A(x) with 
given function F(x) 

A(A(x)) = F(x) 

appeared in 1820 (Babbage, with F(x) = x). 

In this section we consider the case, when the expression for coefficients of a generating 
function A(x) = ^2 n>0 a(n)x n is not known, but its composition with itself is known. 

Theorem 3. Suppose F(x) = ^2 n>0 f(n)x n is a generating function, where /(l) ^ 0, 
F A (n,k) is the composita of F(x), A(x) = Yl n >o a ( n ) xn ^ s ^ e generating function, which 
is obtained from the functional equation A(A(x)) = F(x). Then for the composita of the 
generating function A(x) we have the following recurrent formula 

/(l)t, ^ n = k- 

1 ^F A (n,k)- £ A A (n,m)A A (m,k)] , n > k. (5) 



/(!)'+/(!)' V m=k+l 

Proof. Using Theorem 2, for the case n = k we get 

A A (n,n)A A (n,n) = F A (n,n). 
According to the formula (3), we obtain 

A A (n,n) = f(l)^. 
Using Theorem 2, for the case n > k we get 

n 

AA ( n > m)A A {m, k) = F A {n, k) 

m=k 

or 

n-l 

A A {n,k)A A {k,k) + A A {n,m)A A {m,k) + A A {n,n)A A {n,k) = F A {n,k). 

m=k+l 

Therefore, we have 

n-l 

A A {n,k) (A A {k,k) + A A {n,n)) = F A {n,k) - ^ A A {n,m)A A {m,k). 

m=k+l 

Since A A (n,n) = a(l) n = /(l)^, we obtain the desired formula 



A A (n,k) 



/(l)*, n = k; 

' n-l \ 

F A (n,k)- A A (n,m)A A (m,k) j , n > k. 

m=k+l J 



/(l) 5 +/(1) ^ 



Next we show that the formula (5) is computable. For the case n = 1 the formula is 
computable 

a a (i,i) = y/m. 

By induction, we put that there exist A A (k, I) for k,l ^ re. 

Then we prove that A A {n + 1, m) is computable for m ^ n + 1. 

1. First case m = n + 1. From Theorem 2, we have 

A A (n + l,n+ l)A A (re + l,re + 1) = F A (re + l,n+ 1). 
Then A A (n + 1, n + 1) is computable. 

2. Next case m = n. From Theorem 2, we have 

A A (n + 1, n)A A (n, re) + A A (n + 1, re + l)A A (n + 1, re) = F A (n + 1, re). 
Since ^4 A (re, re) and ^4 A (re + 1, n + 1) are known, A A {n + 1, re) is computable. 

3. The case m = re — 1. From Theorem 2, we have 

A A (re + 1, re - l)A A (re + 1, re + 1) + A A (n, n - l)A A (n + 1, re) + 

+A A (n - l,n - l)A A (re + 1, re - 1) = F A (n + 1, re). 

Since yl A (re + l,re + l),A A (re + l,re) are obtained previously, and other are known by 
the condition, A A {n + 1, n — 1) is computable. 

Therefore, by induction we can obtain all values of A A {n + 1, m) for m ^ n + 1. 

□ 

For applications of Theorem 3 we give some examples. 

Example 4. Let us consider the iterative functional equation A(A(x)) = sin(x) (for details, 
see A048602 , A048603 in [14]). According to [10], the composita of sin(x) is equal to 



2 fc re! 

m=0 

Using Theorem 3, we obtain following recurrent formula 

1, re = k; 

A ^ k)= i ifl±MpX:^)(27re-fcr(-l)^--- £ A A (re,m)A A (m^)V n > k. 

\ m=0 m=k+l J 

Therefore, using the formula (2), the expression for coefficients of A(x) = En>o a ( n ) x " 

is 

a(re) = A A (re,l). 



Example 5. Let us consider the iterative functional equation A(A(x)) = e x — 1 (for details, 
see A052104 , A052122 , and A052105 in [14]). Stanley studied the problem of finding an 
expression for the coefficients of A(x), however the problem remains unsolved [15, ex. 
5.52.c]. 

Next, we obtain the composita of the exponential generating function [10] 

F( x ) = e x - 1. 

Raising this generating function to the power of k and applying the binomial theorem, 
we obtain 



* - if = £ (*V (-i: 



k-j 



Since 



we get 



x\k ^ ^ K 



n 

n>0 



r,!^' 



j=0 

or since the general formula for the Stirling numbers of the second kind is given as follows: 
we have 

f>, (6) 

Here {^} stand for the Stirling numbers of the second kind(see [11, 16]). The Stirling 
numbers of the second kind {^} = S(n, k) count the number of ways to partition a set of n 
elements into k nonempty subsets. 

Using Theorem 3, we obtain following recurrent formula 

{1, n = m, 

ifS {ft" E AHn,m)A*(m,k)), n > k. 
\ m=k+l J 

Therefore, using the formula (2), the expression for coefficients of A(x) = J2 n >o a ( n )^T 

is 

a(n) = n\A A (n, 1). 

Other examples of solutions of iterative functional equations written in the On-line En- 
cyclopedia of Integer Sequences (see Table 1). 



5 



Table 1: Iterative functional equations and their solutions in OEIS 



N of sequence 


Equation 


A030274 


A(A(x)) - 


A097088 


A(A(x)) = x + x 2 


A097090 


A(A(x)) = x(l + 2xf 


A048607 


A(A(x)) = ln(l + x) 


A072350 


A(A(x)) = tan(x) 


A199822 


A{A(x)) - !- 4 --vi-8x- 


A199823 


MA(x)) = ^3 


A199852 


A(A(x)) = xexp(2x) 


A199917 


A(A(x)) = i(2-2cos(z)) 



4 Integer properties of an iterative functional equation 

A(A(x)) = F(x] 



In this section we study an iterative functional equation A(A(x)) = F(x), where F(x) 
Xln>o f( n ) xH is a generating function with integer coefficients. 
First for case /(l) = 1, the formula (5) has form 



n = k 



A A (n,k)-{ if F A^ k ^_ j2 A A (n 1 m)A A (m,k)] , n > k. ^ 

V m=fc+l / 

According to the above formula (7), we present the composita of A(x) in tabular form 



^2.1 



V 2 A i^-<i) j ( 8 ) 

((P2 A l^ A 2-2^ A l)^4 A 3-2^2 A in A 2+8^4 A l) j^Vj^) £^ -| 

16 8 2 

Considering Triangle 8, we can see that not all expressions 2 n ~ k A A (n, k) is integer (for 
example for A A (3, 1)). Now we give the following theorem. 

Theorem 6. Suppose F(x) = ^2 n>0 f(n)x n is a generating function with integer coefficients 
such that /(l) = 1, F A (n, k) is the composita of F(x), A(x) = Yl n >o a ( n ) xn ^ s ^ e generating 
function, which is obtained from the functional equation A(A(x)) = F(x). Then for the 
composita of the generating function A(x) the expression 

A n ~ k A A (n,k) 

is integer. 
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Proof. To prove the theorem we use the formula (7). 

For the case n < 5 the expression A n ~ k A A (n, k) is integer for all k ^ n (see Triangle 8). 
By induction, we put that A n ~ k A A (k, I) is integer for k,l ^ n. 
Then we prove that 

A n+l ~ k A A {n + l,k). 

is integer for all k ^ n + 1. 

According to the formula (7), we have 

A n+1 ~ k A A (n + l,k) = 
= 4 n+i-fcl lp A (n + l,k) - A A (n + l,m)A A (m,k)) = 

\ m=k+l / 

= I j A n+1 ~ k F A (n + 1, k) - 4 n+1 ""M A (n + 1, m)4 m - k A A (m, k) ] . (9) 

\ m=k+l J 

Since 4 n+1_m F A (n + l,m) is even and 4 k ~ m A A (k, m) is integer, we need to prove that 
the expression 

A n+1 - m A A (n + l,m) 

is even for m < n + 1. 

1. First case m = n. According to the formula (9), we have 

AA A (n + 1, n) = AF A (n + 1, n)-. 

Then 4 n+1 ~ m A A (n + 1, m) is even for m = n. 

2. Next case m — n — 1. According to the formula (9), we have 

4 2 A A (n + 1, n - 1) = - (4 2 F A (n + 1, ra - 1) - 4A A (n, n - l)4A A (n + 1, n)) . 

Since above conditions, the expression A n+1 ~ m A A {n + 1, m) is even for m = n — 1. 

3. The case m = n — 2. According to the formula (9), we have 

A z A A (n + l,n - 2) = - (4 3 F A (n + 1, n - 2) - 4 2 A A (n, n - 2)4 1 yl A (n + 1, n)+ 

+A 2 A A (n, n - 2)A l A A (n + l,n) + A 1 A A (n -l,n- 2)4 2 y4 A (n + 1, n - 1)) . 

Since 4 2 v4 A (n + 1, n — 1) is even, and other are known by the condition, the expression 
A n+1 ~ m A A (n + 1, m) is even for m = n — 2. 
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Therefore, by induction we can obtain all values of A A {n + 1, m) for m < n + 1 which 
will be even. 

The theorem is proved. □ 

Corollary 7. Suppose F(x) = ^2 n>0 f(n)x n is a generating function with integer coefficients 
such that f(l) = 1, F A (n, k) is the composita of F(x). Then the generating function A(x) = 
X]„>o a ( n )% n , which is obtained from the functional equation A(A(x)) = , has integer 
coefficients. 

Proof. According to [10], the composita of F ^> is equal to 

4 n - k F A (n,k). 

Next we consider the formula (7) 

1, n = k; 



A A (n, k) 



n-l 

l(4 n - k F A (n,k)- Yl A A (n,m)A A (m,k) ) , n > k. 

m=k+l 



Replacing A A (n, k) by the composita 4" k B A (n, k), we get 



A n - k B A (n,k) = - (A n - k F A (n,k) - ^ m B A {n,m)A m - k B A {m 1 k) \ . 

\ m =k+l J 

According to Theorem 6, the expression A n - k B A (n,k) is inte ger. 
Therefore, the generating function A(x) such that satisfies to the equation 

A{A{x)) = 

where F{x) = ^2 n>0 f(n)x n is a generating function with integer coefficients and /(l) = 1, 
has integer coefficients. 

The corollary is proved. □ 

Example 8. Let us consider the iterative functional equation 



A(A(x)) = y — . 

Also the generating function A(x) satisfies the following functional equation (for details, see 
A213422 in [14]) 

A(A(x) - 4A{x) 2 ) = x, A{A{x)) = x + 4A{A{x)) 2 . 
Here F(x) is the generating function for the Catalan numbers 



l-y/l-4x 

b (x) = . 



According to [10], the composita of F ^ = 1 vl 16 x . j s equal to 

A n ~ k F A (n, k) = 4 n - fc ( 2n ~^~ 

Using Theorem 3, we obtain following recurrent formula 

1, n = k; 

A/„ ;„\ _ ) / Ti-l 



V fe=m-l / 



According to the above formula, we present the composita of v4(x) in tabular form 
1, 

2F^i, 1 
8F A 3,1-2F 2 A A 3>2 , 2F 3 A 2 , 1 

1^^2,1^3,2 or 3, 1/^4,3 °- r 2, 1^4,2 ^ oz - r 4,H °- r 4,2 z - r 3,2 r 4,3' zr 4,3' 1 

Therefore, using the formula (2), the expression for coefficients of A(x) = J2 n>0 a(n)x n is 

a (n) = A A (n,l). 

The first elements of the obtained generating function A(x) are shown below 
x + 2x 2 + 12a; 3 + 96a; 4 + 880a; 5 + 8720a; 6 + 90752a; 7 + • • • 

5 Iterative functional equation A 2n {x) = F(x) 

We start with the case A(A(A(A(x)))) = F(x). Let us get the following iterative functional 
equation 

B(B(x)) = F(x), 

where B(x) = A(A(x)). 

Using Theorem 3, we obtain the composita of B(x) 



B A (n, k) 



/(l)*, n = k; 

» (F A (n,k)- E B A {n,m)B A {m,k)Y n > k. 

/(1)TT+/(1)2 y m =k+l J 



Then using the formula (5) again for the equation B(x) = A(A(x)), we obtain the 
composita A A (n, k). 

Let us consider this approach in the following example. 

Example 9. Suppose we have the following iterative functional equation 

A(A(A(A(x)))) =x + 16x 2 . 
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According to [10], the composita of F(x) = x + 16x 2 is equal to 

F A (n,k) = ( k )l6 n ~ k . 



y n — k / 

For the generating function B{x) such that B(B(x)) = x + 16x 2 the composita is equal to 

1, n = k; 

, n > k. 



B A (n,k) 



n-1 
i=fe+l 



Since B(x) = A(A(x)), we obtain the composita of A(x) 

1, n — k; 



A A (n,k) 



n-1 



B A (n,k) - £ A A (i,k) A A (n,i) 



, n > k. 



i=k+l 

Therefore, using the formula (2), the expression for coefficients of A(x) = J2 n >o a ( n ) xn ^ s 

a(n) = A A (n, 1). 

The first elements of the obtained generating function A(x) are shown below 

x + Ax 2 - 48x 3 + 960x 4 - 23296x 5 + 616448x 6 - 16830464a; 7 + • • • 
For details, see A141119 in [14]. 

Generalize the above proposed approach for the iterative functional equation A 2n (x) = 
F(x). 

1. Make substitution B(x) = A^ 2 " \x) and solve B(B(x)) = F(x). Then we obtain 
B A (n,k). 

2. Next make substitution C(x) = B^ 2 \x) and solve C{C{x)) = B{x) with respect to 
C(x). Then we obtain C A (n, k). 

3. Making the analogic substitution until n — 1, we obtain the desired solution of A(x). 
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